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SUMMART 


This  report  Is  concerned  with  a  linear  time  varying  approximation  to  the 
dynamics  of  low-speed  flying  machines.  Simplifications  and  approximations 
are  applied  widely  in  order  to  emphasize  essential  aspects . 

The  range  of  time  v6urlation  is  described  in  terms  of  frozen  system.' lodi  of 
the  roots  corresponding  to  the  predominant  mode  of  a  system.  The  rate  of 
the  time  variation  is  described  In  terms  of  the  deviation  from  the  frozen 
system  approximation.  An  analog  computer  study  was  made  to  specify  qmnti- 
tatlvely  those  rates  of  time  variation  which  cannot  be  considered  as  slow. 

The  longltiodinal  dynamics  of  VTOL  aircraft  is  stvidled  as  an  exan5>le  in 
rather  general  terms.  Approximations  and  the  application  of  root  locus 
methods  in  terms  of  the  most  significant  stability  derivatives  lead  to  a 
construction  describing  the  behavior  of  the  oscillatory  roots  during 
transition. 

The  results  are  rased  in  a  discussion  of  the  following  variable  feedback 
configrorations;  direct  feedback  adjustments,  adaptive  feedback,  and  pro¬ 
grammed  feedback  adjustments. 
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CQNCHJSIOWS 


Two  distinct  aspects  of  time  varying  systems  are  the  range  and  the  rate  of 
time  variation.  The  q.uestion  of  when  the  time  variation  shoiild  be  considered 
"fast"  is  answered  on  the  basis  of  practical  considerations  by  means  of  a 
computer  study.  The  rate  of  time  variation  should  be  considered  fast  if, 
during  the  time  of  one  period  of  the  predominant  mode,  the  spring  constant 
variation  is  eqml  to  (or  more  than)  and/or  the  damping  coefficient  vari¬ 
ation  is  equal  to  (or  more  than)  .5*  In  these  cases  the  system  response  is 
markedly  different  from  that  obtained  by  the  frozen  system  approximation. 

The  range  of  the  time  variation  and  the  behavior  of  the  predominant  roots  can 
be  predicted  to  a  large  extent  by  a  construction  based  on  simplifying  as¬ 
sumptions.  This  method  makes  the  tracking  of  the  influence  of  important  sta¬ 
bility  derivatives  very  straightforward.  With  this  method,  considerable 
a  priori  knowledge  can  be  provided  on  the  behavior  of  the  oscillatory  roots 
during  VTOL  transition.  This  knowledge  is  vised  in  a  dlscvission  of  three 
different  principles  which  can  be  used  for  artificial  stability  augmentation 
by  means  of  feedback.  The  advantages  of  programmed  feedback  gain  adjustment 
versvLS  adaptive  control  are  enqjhaslzed. 

Recommendations . 


The  performance  of  an  artificial  stability  augmentation  system  with  programmed 
gain  adjustments  depends  largely  on  the  amount  of  information  available. 
Therefore,  further  refinements  of  the  presented  approximating  methods,  ea- 
pecially  in  mid-transition,  vould  be  of  great  value.  The  full  potential  and 
limitations  of  such  a  system  shovild  be  investigated  by  exploring  the  sensi¬ 
tivity  of  the  performance  to  deviations  from  assumed  design  parameters. 
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mTHODUCnCK 


One  common  characteristic  of  low-speed  flying  machines  is  their  ability  to 
achieve  rapid  changes  of  the  flight  condition.  The  accent  is  on  the  word 
"rapid,"  and  the  present  report  is  concerned  primarily  with  this  aspect  of 
the  flying  characteristics  of  such  vehicles . 

The  reason  why  low-speed  flying  machines  exhibit  this  characteristic  can  be 
made  plausible  by  considering  Newton's  law:  P  =  ma.  If  a  certain  accelerat 
Ing  force,  F,  is  acting  on  a  body  for  a  time  interval,  t,  the  achieved  ve¬ 
locity  can  be  expressed  in  this  form: 


The  smaller  the  initial  velocity,  V^,  the  larger  is  the  relative  change  in 
velocity.  Most  of  the  aerodynamic  stability  and  control  derivatives  vary 
with  therefore,  a  set  of  linearized  differential  eq.uatlons  can  be  ex¬ 
pected  to  become  invalid  more  rapidly  at  low  speeds  than  at  high  si)eeds. 

It  should  be  noted  here  that  the  rapid  changes  of  characteristics  at  very 
high  speeds,  although  similar  in  effect,  are  of  a  different  origin.  The 
very  high  velocities  enable  the  vehicle  to  encounter  rapid  changes  of  the 
environment;  for  example,  of  the  air  density. 

In  both  cases,  at  very  low  and  at  very  high  velocities,  it  is  actually  the 
nonlinear  variation  of  the  dynamic  pressure  which  causes  the  conventional 
set  of  linear  eq,uatlons  to  become  invalid,  therefore,  the  findings  of  this 
report,  although  presented  with  low-speed  flying  machines  in  mind,  are  appli 
cable  to  dynamic  problems  at  both  extremes  of  the  velocity  range. 

A  rigorous  approach  to  the  outlined  problem  would  lead  to  a  set  of  nonlinear 
differential  eq.viations  valid  over  the  required  velocity  range.  In  the 
present  state  of  the  art  such  an  approach  could  not  be  expected  to  yield 
useful  results.  A  more  promising  approach  tries  to  extend  the  classical 
method  of  small  perturbations  so  that  approximating  results  co\ild  be  ex¬ 
pected  to  reveal  the  essential  behavior  of  the  nonlinear  system. 

Consider  a  general  physical  system  as  described  by  a  set  of  differential 
equations.  There  are  three  kinds  of  "elements"  which  can  be  distinguished 
in  most  differential  equations: 
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1.  Functions  of  the  Independent  variables  onJy  (inputs). 

2.  Dependent  variables  and  their  derivatives  (outputs). 

3.  Coefficients  (parameters). 

The  equivalent  terms  generally  \jsed  are  in  parentheses .  A  simplified  repre¬ 
sentation  of  a  fairly  general  system  can  be  given  in  the  foUovlng  form; 


This  somewhat  unconventloneil  representation  has  been  chosen  in  order  to  empha¬ 
size  the  general  nature  of  the  relationships.  A  nonlinear  differential 
equation  can  be  considered  a  variation  of  an  approximating  linear  differential 
equation  in  which  the  coefficients  are  ftmctlons  of  the  dependent  variables. 
The  line  connecting  the  outputs  with  the  parameters  then  represents  the  es¬ 
sential  characteristic  of  nonlinear  systems  which  is  the  dependence  of  coef¬ 
ficients  on  the  dependent  variables.  Linearization  means  cutting  this  line  as 
indicated,  l.e.,  considering  the  coefficients  to  be  Independent  of  the  de¬ 
pendent  variables.  Total  remcrval  of  all  lines  affecting  the  parameters  is 
equivalent  to  a  further  simplification  to  a  constant  parameter  system  de¬ 
scribed  by  differential  equations  with  constant  coefficients.  This  latter 
approach  is  \ased  when  small  perturbations  of  an  equlllbrlxmi  (trim  condition) 
are  investigated.  In  this  case  there  is  a  conventional  llneeu*  system  with  ot 
withovit  feedback. 

The  effect  of  the  loop.  Indicated  by  the  line  between  outputs  emd  parameters, 
on  such  a  linearized  system,  will  now  be  discussed.  From  the  viewpoint  of 
dynamic  response#  there  are  two  distinct  effects  of  the  nonlinear  relation¬ 
ships  represented  by  the  upper  loop  in  the  sketch. 


1.  Nonlinear  effects  may  chemge  the  dynamics  in  the  vicinity  of  the 
eoLviilibrivnn  so  that  linearized  perturbation  equations  are  not  valid.  For 
example,  satiiration,  backlash,  etc.,  cavise  such  effects.  Several  methods 
(e.g.,  phase  space,  describing  f met ions)  can  be  used  to  deal  with  these 
problems . 

2.  Nonlinear  effects  may  change  the  eq.viil1brium  (trim)  conditions  to 
such  an  extent  that  the  coefficients  of  the  linear  pertvurbatipn  equation 
must  be  changed.  For  example,  the  change  of  dynamic  pressure  or  the  change 
of  mass  with  fuel  consumption  came  such  effects.  It  is  this  type  of  ef¬ 
fect  with  which  the  present  report  is  concerned. 

It  must  be  emphasized  that  the  above  distinction  is  not  based  on  mathemati¬ 
cal  but  rather  on  practical  considerations.  Actually,  the  same  non¬ 
linearity  may  give  rise  to  both  types  of  effects,  depending  upon  the  forces 
acting  on  the  system.  The  first  case  ocems  if  the  dynamic  motion  result¬ 
ing  from  distiurbances  and  control  actions  involves  the  nonlinearity;  however, 
the  equations  describing  the  dynamics  remain  unaltered  throughout  the  regime . 
The  problems  arising  from  this  situation  are  not  of  concern  here.  It  is 
assumed  throvighout  the  present  report  that  the  "small  perturbation  equations" 
of  the  system  are  linear. 

The  second  case  occurs  when  an  average  dlsttirbance  or  control  action  moves 
the  equilibrium  state  along,  the  nonlinearity.  The  change  of  the  linear 
term  in  the  Taylor  series  expansion  of  the  nonlinearity  means  continuous 
changes  in  the  coefficients  of  the  small  pertvurbatlon  equations.  This  con¬ 
sideration  leads  to  the  following  linearization  of  the  nonlinear  problem. 

It  is  assumed  that  the  average  forces  and  moments  are  continuously  balanced 
(trimmed) .  Small  deviations  from  the  trim  condition  appear  as  forcing  terms 
in  the  perturbation  equations.  If  thd  assumption  is  now  made  that  the  time 
history  of  the  trim  conditions  can  be  predicted,  then  the  parameters  of  the 
pert\irbatlon  equations  can  be  considered  functions  of  time  instead  of 
fimctlons  of  the  system  output.  With  these  asstmptions  the  nonlinear  system 
is  approximated  by  a  linear  time  varying  system.  The  upper  loop  in  the 
Figure  on  page  4  has  been  cut  and  the  arrow  pointing  to  the  parameters  repre¬ 
sents  a  function  of  the  time,  t,  which  can  be  considered, in  a  generalized 
way,  as  an  additional  parametric  input  to  the  system. 

If  the  predicted  time  variation  is  slow  so  that  the  parameters  of  the  per- 
tvirbatlon  equatlcm  can  be  considered  constant  for  the  duration  of  at  least 
one  transient  response,  then  the  system  can  be  further  sinqpllfied  to  a 
"quasl-stationary"  system.  In  this  case,  loci  of  the  poles  and  zeros  can 
be  plotted  with  time  as  the  rtmnlng  parameter.  Such  a  plot  represents  the 


5 


so-called  "frozen  system"  approximation  to  the  time  variable  system.  The 
somewhat  vague  fashion  in  which  the  "quasl-statlonary"  or  "frozen"  concept 
is  introduced  indicates  the  difficulty  in  defining  a  border  line  be+-ween 
time  variable  and  constant  parameter  systems.  Jigain,  the  distinct iv>n  has 
to  be  based  on  practical  engineering  rather  than  on  mathematical  considera¬ 
tions  .  Obviously,  problems  arise  when  the  time  variation  is  such  that  the 
quasi-statlonary  approximation  does  not  hold. 

This  report  is  concerned  with  the  linear  time  variable  approximation  to  the 
dynamic  problems  of  low -speed  flying  machines.  Only  smooth  and  monotonlc 
time  variation  axe  considered.  Simplifications  and  approximations  are  ap¬ 
plied  widely  in  order  to  emphasize  essential  aspects  rather  than  detailed 
calculations.  .One  major  aspect  investigated  is  the  effect  of  the  speed  of 
parameter  variation.  Various  feedback  stability  augmentation  configurations  - 
are  considered.  As  an  example,  the  transition  of  a  VTOL  aircraft  is  dis- 
ciissed  in  more  detail. 
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THE  FRCEEN  SYSTEM  AEPRODOMATION 


Time  varying  linear  systems  can  be  described  by  linear  differential 
equations  in  which  the  coefficients  are  functions  of  time.  If  the  time 
variation  is  slow,  the  coefficients  will  not  vary  to  any  great  extent  with¬ 
in  a  given  time  interval.  The  definition  of  "slowness"  involves  the  speci¬ 
fication  of  the  time  interval  to  be  considered  sis  well  as  the  specification 
of  the  coefficient  variation  to  be  tolerated.  For  practical  reanons  these 
specifications  have  to  be  determined  with  the  act\jal  output  of  the  system 
in  mind. 

A  stationary  system  is  usually  described  either  in  the  frequency  domain 
(transfer  fimctlon  or  amplitude  and  phase  characteristic)  or  in  the  time 
domain  (weighting  fvinction  or  impulse  response).  The  frequency  domain  ap¬ 
pears  to  be  the  less  natural  choice  in  which  speed  of  time  variation  should 
be  defined  since  it  implies,  by  definition,  steady  state  test  frequency 
signals.  The  impulse  response,  on  the  other  hand,  provides  the  significant 
information  within  a  limited  time  interval.  A  time  Interval,  at  the  end  of 
which  it  is  possible  to  make  a  judgement  on  the  essential  characteristics 
of  the  system,  seems  xoseful  as  a  basis  for  the  investigation  of  the  rate  of 
time  variation.  The  significance  of  the  time  variation  can  be  J\adged  by  the 
deviation  from  the  response  of  a  stationary  system.  These  ideas  are  now 
considered  in  more  detail. 

Assume  a  dlfferentlail  equation  of  a  time  varying  system  with  an  Impulse  ap¬ 
plied  at  the  time,  T,  in  the  followicg  form: 

2^  fj^(t)  ^  X  (t)  =  6(t)  (1) 

i  =  o 

If  the  variation  of  the  coefficients,  fjj^(t),  is  very  slow,  then  the  quasl- 
stationary  equation 

1  =  0 

is  a  good  approximation.  If  the  roots  of  the  corresponding  characteristic 
equation  are  plotted  with..  T  as  the  running  parameter,  the  frozen  system 
loci  of  the  time  varying  system  are  obtained.  In  general,  the  time  varia¬ 
tion  of  the  zeros,  as  veil  as  of  the  poles  of  a  transfer  function,  can  be 
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characterized  in  this  way.  The  loci  do  not  change  if,  t.  Is  replaced  by,  tjt. 
The  time  variation  is  faster  when  Increases.  Obviously,  the  approxi¬ 
mation  becomes  worse  with  increasing,  tj.  The  frozen  shlution,  y(t),  can  be 
used  as  a  reference  with  which  the  solution,  x(t),  of  the  time  variable  can 
be  compared. 

The  problem  is  to  define  a  time  interval  which  permits  a  meaningful  charac¬ 
terization  of  the  frozen  system  solution  eind  which  can  also  be  used  for 
evaluating  the  deviation  of,.  x(t),  from  the  frozen  system  solution;  y(t). 
Solving  this  problem  in  its  full  generality  will  not  be  attempted.  Use  will 
be  made  of  two  restrictions  determined  by  the  subject  covered. in  this  report; 
(l)  to  consider  systems  in  which  the  time  variation  is  smooth  and  mono¬ 
tonic,  at  least  piecewise;  (2)  to  examine  the  effect  of  the  time  variation 
from  the  point  of  view  of  its  effect  on  a  feedback  control  applied  to  the 
time  varying  system.  The  first  restriction  permits  some  highly  simplifying 
approximations;  the  second  restriction  is  of  great  help  in  determining  a 
time  interval  for  a  meaningful  characterization  of  the  rate  of  time  vari¬ 
ation  and  its  effects. 

A  discussion  of  some  consequences  of  feedback  control  is  presented  first. 

The  frozen  system  is  considered  then  the  deviations  due  to  fast  time  vari¬ 
ation.  As  mentioned  before,  the  frozen  system  can  be  represented  by  loci 
of  poles  and  zeros  on  the  complex  plane.  At  each  instant,  t,  determines  a 
certain  pole-zero  configuration  or  pulse  response.  A  pulse  response  of  the 
frozen  constant  parameter  system  is  the  sum  of  component  response  modes. 

The  modes  are  determined  by  the  poles;  the  relative  amplit\ides  and  phases 
are  determined  by  the  relative  location  of  the  zeros  and  poles .  Consider¬ 
ing  the  resjxDnse  as  the  input  to  a  feedback  branch  (automatic  or  human),  it 
is  important  to  distinguish  between  dominant  modes  and  modes  of  minor  sig¬ 
nificance  .  Even  resjKsnses  of  systems  of  high  order  can  usually  be  matched 
qviite  will  with  a  second,  third,  or  possibly  fourth-order  system  response. 
Th\is  in  the  case  of  a  time  variable  system  the  variation  of  the  dominant 
modes  is  of  primary  significance. 

It  is  very  dlfflc\ilt  to  be  quantitative  about  a  general  criterion  for  de¬ 
termining  whether  or  not  a  mode  should  be  considered  predominant.  However, 
this  study  is  restricted  to  aircraft  configurations.  Equations  of  higher 
than  fo\irth-order  need  not  be  considered.  The  relative  importance  of 
different  modes  should  be  Judged  on  the  basis  of  relative  asplltudes. 

The  relative  amplitudes  of  the  different  modes  in  the  pulse  response  can  be 
determined  graphically  frcmi  the  pole-zero  plot  of  the  transfer  function 
(Reference  l).  The  ampllt\ide  for  any  individual  pole  is  determined  by 
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dividing  the  product  of  the  distances  of  the  pole  from  the  zeros  by  the  pro¬ 
duct  of  the  distajcces  from  the  other  poles .  The  amplltiode  of  an  oscilla¬ 
tory  mode,  represented  by  a  pair  of  conjxigate  complex  poles,  is  twice  the 
amplitude  determined  for  one  of  these  poles.  Any  arbitrary  scale  can  be 
chosen  since  only  the  relative  amplitudes  are  of  interest. 

The  following  sketch  shows  some  significant  points  in  this  procedvire.  Con¬ 
sider  the  following  pole-zero  configuration. 


When  the  amplitudes  of  the  oscillatory  mode  and  the,  mode  are  comgpared 
in  the  form  of  a  ratio,  the  distance  between.,  p^.  (or  pij.)  and,  p^,  drops  out. 
The  ratio  can  be  written  as 
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In  this  example,  Is  at  a  similar  distance  from  both  poles.  The  ratio  of 
the  distances  to,  p^,  together  vith  cos  e^,  result  in  a  factor  of  approxi¬ 
mately  two  in  favor  of  The  decisive  factor  is  caused  by  the  fact 

that  Zj^  is  much  closer  to^p^  than  to  Po.  Therefore,  it  can  be  safely  esti¬ 
mated  that  the  amplitude  of  the  oscillatory  mode  is  more  than  five  times  the 
amplitude  of  the  p^  mode . 

The  relative  amplitude  ratio  A2/A2  of  the  P2_  and  P2  modes  is  determined  in 
this  example  primarily  by  the  square  of  the  ratio  of  the  distances  to  p^ 

(or  pi^)  because  the  arrangement  of  the  zeros  is  fairly  symmetrical  with  re¬ 
spect  to  the  two  real  poles .  Therefore 

I  '  2 

,  A2 I  cos  £2 

'  -^li  cos^ 


(M 


The  considerations  in  this  example  can  be  utilized  generally.  For  example, 
the  conditions  imder  which  the  short  period  mode  of  a  pitch  response  is  pre¬ 
dominant  can  be  easily  estimated.  Since  no  objective  criterion  for 
"predominant"  can  be  established,  some  r\ile  of  thumb  must  be  used.  The 
amplitude  of  a  predominant  mode  should  not  be  less  than  four  to  five  times 
the  amplitude  of  any  other  mode.  This  is  a  reasonable  criterion  for  a 
fourth  order  system. 

NORMALIZATION 


It  is  convenient  to  use  a  normalization  so  that  the  effect  of  time  variation 
on  a  particular  mode  can  be  investigated  in  a  general  form.  In  order  to  in¬ 
vestigate  some  fundamental  aspects  of  the  time  variation,  only  one  pair  of 
complex  conjugate  roots  is  considered.  This  is  equivalent  to  replacing, 
for  example,  an  aircraft  pitch  attitude  response  with  the  short  period  oscil¬ 
latory  mode.  This  oversimplification  will  be  helpful  in  the  investigation 
of  practical  effects  of  fast  time  variation.  The  effect  of  deviations  from 
this  simplified  model  is  dlscxissed  later. 

The  homogeneous  differential  eqiiation  of  the  simplified  time  varying  system 
of  second  order  is; 


(5) 

The  solution  in  the  vicinity  of  a  frozen  time  solution  is  now  examined. 
Using  a  conventional  notation  of  fixed  parameter  systems 


x  +  f2^(t)x  +  f2(t)x  =  0 
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■where  the  ■’ndex 


V’)  - 

F;.  stands  for  "frozen.” 


In  the  course  of  transitions,  the  time  variations  of  parameters  are  not  only 
smooth  hut  also  monotonic  over  the  full  range,  or  at  least  over  large  por¬ 
tions  of  the  range.  Therefore,  it  is  reasonable  to  investigate  the  effect 
of  time  variation  by  considering  only  the  first  tvo  terms  of  the  Taylor 
series  of  the  coefficient  time  functions:  the  time  variable  damping  and 
spring  coefficients. 


r. 


+  I^^F^Fo  +  ■  ■'^>1  ^  + 


+  D(t  -  t) 


X  =  6(t) 


(6) 


For  convenience,  the  origin  of  time  js  shifted  to  t.  Using  the  notations 
B  =  p(4,  D  =  wp  =  \l  1  - 


X  +  +  (a^^+  iia^t)x  =  5(o) 


(7) 


Normalization  is  obtained  by  a  linear  transformation 

Cf^o^ 

y  =::xe 

and  a  change  in  time  scale 


f  =  Wpt 
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with 


P  =  P 

Q  =  |jl  +  p  tan  \|r  (8) 

the  normalized  eq.uation  is 

y  +  Pt'y  +  (l  +  ftt')v  =  &(o)  (9) 

Here  t'  is  dimensionless  time  as  defined  above.  This  normalization  puts 
the  frozen  poles  into  the  (+J,  -j)  points  of  the  Sy  plane.  The  effect  of 
time  variation  can  be  investigated  in  this  plane.  The  simple' transfor¬ 
mation  relationships  are  described  in  more  detail  in  Appendix  I. 

The  solution  of  the  differential  eq.uation  with  linearly  time  varying  coef¬ 
ficients  can  be  obtained  in  the  form  of  hyijergeometric  functions.  However, 
for  engineering  purposes,  such  a  solution  in  itself  is  not  very  helpful  in 
the  evaluation  of  fast  time  variation.  An  analog  computer  stu^  was  pre¬ 
ferred  beca\ise  of  the  extreme  flexibility  offered  by  this  approach.  High 
accuracy  is  not  Important  because  of  the  approximations  Involved  in  the 
study.  The  goal  is  rather  to  identify  significant  effects  of  time 
variation. 
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ANALOG  COMPUTER  STUDY 


The  scheme  for  experiments  with  the  normalized  second  order  system  was  deter¬ 
mined  as  follows.  Since  the  normalization  trajisforms  straight  lines  into 
straight  lines.  It  seemed  advantageous  to  explore  time  variations  which  move 
the  roots  along  straight  lines,  rather  than  purely  linear  time  variations 
which  move  the  poles  along  circles. 


The  normalized  system  equation 


+  (1  •  Ci^t'  +  Q2t'^)y  =  0 

(10) 

for  "frozen"  roots; 

Pt'  \ 

p2t 

(11) 

1  +  ^2^*^  *  +  ^2*^ 

-  “TJ— 

If  Q.2-.  - 


+  Q,f 


,  this  locus  becomes  a  pair  of  straight  lines; 


S^t‘)  =  -  ^  t  J(1  + 

-1 

The  slope  tan  and  the  constant  speed  of  variation 


dt*^ 


+  Ql 


(12) 


are  the  same  as  those  determined  for  the  linearly  time  varying  case  at 
t'  =  0  in  Appendix  I.  It  was  ^decided  that  various  speeds  of  variation  along 
eight  different  directions  45  from  each  other,  originating  in  the  "frozen" 
pole  pair  (+J,  -j)  should  he  investigated. 


It  was  rather  difficult  to  choose  a  basis  for  the  evaluation  of  the  effect 
of  time  variation.  In  the  case  of  slow  time  variation  the  frozen  system  is 
a  good  approximation.  The  goal  of  the  analog  coji^Titer  study  was  to  explore 
the  deviations  from  the  frozen  system  solution  due  to  fast  time  variation, 
from  the  viewpoint  of  the  application  of  feedback.  As  stated  above,  in 
order  to  specify  a  speed  of  time  variation  a  time  Interval  and  a  criterion 
for  the  comparison  must  be  chosen.  The  concepts  of  apparent  dan^jing  and 
apparent  frequency  change  were  used  in  this  choice.  Although  theoretically 
an  "instantaneous"  dancing  and  frequency  could  be  defined,  these  concepts 
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-have  no  practical  mesining  because  some  finite  time  interval  is  necessary  to 
measure  damping  and  freq.uency.  A  human  pilot  is  not  very  sensitive  to 
slight  frequency  changes;  damping  can  be  Judged  most  easily  by  i)eak  ratios; 
therefore,  at  least  one  half  cycle  should  be  iised  for  comparison.  In  order 
to  provide  a  mechanical  method  for  the  comparison,  the  apparent  damping  and 
frequency  were  assumed  to  be  the  damping  and  frequency  of  an  adjustable  con¬ 
stant  parameter  system.  The  response  of  this  system  with  an  identical  in¬ 
put  was  matched  to  the  response  of  the  time  variable  system.  The  match  was 
made  in  the  "least  Integral  square  error"  (ISE)  sense.  This  procedure  was 
set  up  on  an  analog  computer.  It  was  found  that  using  only  one  half  cycle 
of  the  frozen  system  response  for  comparison  was  not  selective  enough; 
therefore,  one  full  cycle  of  the  undamped  frozen  system  response  was  chosen 
as  a  courparison- Interval. 

A  block  diagram  of  the  coii5>uter  setup  is  shown  in  Figure  1.  The  equation 
of  the  normalized  time  variable  system  is  given  by 

y%  Pty  +  (1  +  =  0  (13) 

where  Q2  chosen  so  that  the  roots  moved  along  straight  loci. 

The  constant  coefficient  system  used  to  approximate  the  time  varying  system 
was 

Va  +  +  (l  +  -  0  (l4) 


where  Kf  and  were  adj\istable  parameters.  Experiments  at  different 

speeds  or  variation  were  conducted  by  applying  identical  Inputs  to  both 
systems  and  measuring  the  displacement  error; 

e(t)  =  y(t)  -  y^Ct) 


This  error  was  squared  by  a  function  generator  and  then  Integrated.  The 
parameters  K.  and  K^j  were  adjusted  to  find  the  partlciilar  comiblnation  which 
minimized  the  ISE  over  one  period  of  the  normalized  frozen  system  response, 
2n  seconds . 


Various  speeds  of  variations  were  used  for  P  and  ^  ,  varying  from  l/kn 
to  2/jt  for  increasing  sprlisg  and/or  damping  coefficifet;  and  from  l/8n  for 
decreasing  spring  and/or  damplz^g  coefficient.  Obviously,  the  greater  the 
speed,  the  larger  is  the  ISE.  Speeds  less  than  1/8k  were  barely  noticed 
(K^  and  K^,  are  less  than  .05).  Speeds  greater  than  l/«  (3OO56  of  spring 
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constant  variation)  generated  huge  errors.  Changing  this  time  Interval 
little  effect  on  Kf  and  for  minimum  ISE  at  slow  speeds .  At  high  speeds 
(1/n,  2/*)  and  ^  were  very  sensitive  to  such  changes.  An  pvantplp'  Indi¬ 
cating  the  selectivity  Is  shown  In  Figure  2. 

Since  very  low  speeds  (P,  l/8jt)  were  found  to  be  approximated  well  by 

the  frozen  system  solution  and  very  high  speeds  can  be  omitted  because  of 
practical  considerations^  the  following  two  speeds  were  selected  for 
comparison: 

P,  Qi  =  l/iht 

P>  Qjl  =  l/2n 

The  speed  =  l/lftt  Increased  the  spring  constant  by  5056  In  2jt  seconds  (the 
period  of  the  frozen  system) .  The  speed  of  =  l/2jt  doubles  the  spring 
constant  In  2n  seconds.  The  speeds  P  =  l/4n  andP  =  l/2x  mean  that.  If 
Ql  =  0,  a  frozen  system  damping  factor  of  .5  and  .707,  respectively.  Is 
reached  after  2n  seconds . 

Flg\ires  3  and  4  show  the  apparent  roots  at  these  speeds  In  different  di¬ 
rections.  The  tips  of  the  vectors  Indicate. the  frozen  system  roots  at  2jc 
seconds.  Figure  5  combines  these  two  figures  by  normalizing  to  the  vector 
lengths . 

These  data,  combined  with  data  obtained- with  other  speeds,  are  plotted  In  a 
different  way  In  Flgiore  6.  The  heavy  line  In  this  flg\are  Indicates  the 
boundary  between  those  values  of  P  and  vhlch  result  In  Ki  0  and  those 
which  result  In  ^  0.  For  the  range  of  speeds  of  condem  It  can  be  seen 
that  the  apparent  "destabilizing"  effect  of  negative  P  can  be  cancelled 
by  a  ixjsltlve  Q  of  twice  the  magnitude.  For  large  values  of  negative  Q, 
the  dynamic  effect  will  be  a  decrease  of  apparent  damping  no  matter  how 
large  P.  (This  Is  meaningful  only  for  Qi  >  -  l/2„  since  If  -  l/2jt, 

the  roots  become  real  within  2jt  seconds ) . 

In  order  to  estimate  the  effect  of  fast  parajneter  variation  for  oscillatory 
roots  emywhere  In  the  complex  plane,  the  Inverse  of  the  normalizing  trans¬ 
formation  must  be  performed.  The  apparent  damping  found  In  the  normalized 
form  estimates  the  change  In  the  damping  of  the  frozen  system  oixfcput.  This 
means  that  the  larger  the  frozen  damping  coefficient,  the  less  the  apparent 
deviation  ca\ised  by  a  fast  time  variation.  Therefore,  It  can  be  stated 
that  the  region  of  Interest  for  this  effect  Is  that  of  low  positive  andnega- 
tlve  damping.  In  those  phases  of  a  VTOL  transition  where  the  oscillatory 
roots  are  In  the  vicinity  of  the  Imaginary  axis,  it  can  be  assxjmed  that 
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deviations  from  the  frozen  system  assumption  will  influence  the  handling 
qiialities  and  the  x)erformance  of  autopilots  designed  on  a  frozen  system 
basis  if  the  parameter  variation  is  fast.  It  has  been  shown  that  fast 
variation  of  the  spring  constant  means  a  to  lOC^  change  within  the 
time  of  one  frozen  system  period.  A  decreasing  spring  constant  causes  more 
deviation  from  the  frozen  system  than  an  increasing  spring  constant.  Fast 
variation  in  damping  means  a  variation  of  .5  bo  .7  in  the  damping  ratio 
during  one  frozen  system  period.  If  a  combination  of  spring  and  damping 
variation  takes  place,  the  effects  may  strengthen  (or  weaken)  each  Other 
with  the  res^Ilt  that  a  smaller  (or  larger)  speed  of  variation  of  the  indi¬ 
vidual  parameters  can  cause  a  similar  effect  on  the  apparent  damping. 

For  VTOL  transitions  the  significance  of  these  findings  is  that  the  effects 
ca\ised  by  fast  parameter  variations  are  expected  to  be  stronger  when  fre¬ 
quency  and  damping  are  decreasing  rapidly  in  the  vicinity  of  the  imaginary 
eixis  than  when  the  roots  move  in  the  opposite  direction  along  a  frozen 
root  locus. 
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AFPRCPaMATIONS  TO  THE 
LOWGrrUDINAL  DYNAMICS  OF  A  VTOL  AIRCRAET 

The  simplified  linearized  rigid-body  eq^uatlons  which  describe  "the  perturbed 
longitudinal  airplane  motion  about  a  straight  and  level  flight  path  may  be 
written  as  follows: 

(S  -  X^)  +  XyW  +  gG  =  0  (16) 

-  z^u  +  (s  -  z^)w  -  vse  =  0 

-  M^u  -  M^,w  +  (S  -  M^)Se  =  M56 

u,  w,  0,  are  the  smaJl  variations  of  the  perturbed  flight  variables  arovind 
the  equilibrium  point.  The  coordinate  system  is  fixed  to  the  body.  The 
X-axis  is  taken  positive  forward  and  is  initially  set  horizontal.  The 
Z-axls  is  positive  downward.  The  only  input  considered  in  these  equations 
is  the  control  moment  around  the  pitch  axis  of  the  airplane. 

1.  MEAR-HOVERING  FUGHT 

The  description  of  the  longitudinal  dynamics  in  near-hovering  flight  can  be 
simplified  by  the  assumption  that  the  vertical  degree  of  freedom  is  un¬ 
coupled  from  the  other  two  (the  horizontal  and  the  angular)  degrees  of  free¬ 
dom.  This  is  Justified  as  long  as  the  stability  derivatives  X^  (variation 
of  horizontal  force  with  vertical  velocity)  and  Mj,  (variation  of  pitch 
moment  with  vertical  velocity)  are  small. 

These  assumptions  are  commonly  used  for  satisfactorily  approximating  the 
modes  of  motion  of  helicopters  near  hovering.  For  VTOL  machines  M^,  can 
be  expected  to  remain  negligible  for  a  larger  velocity  range  near  hovering 
them  is  the  case  for  conventional  single  rotor  helicopters,  because  the  de¬ 
stabilizing  contrib\ition  by  the  rotors  which  increases  as  the  VTOL  airplane 
gains  speed  is  balanced  to  some  extent  by  the  tail  effect.  Moreover,  for 
VTOL  machines  of  the  tilt-wing  type,  the  effect  of  X^,  on  the  treuisfer 
function  is  expected  to  be  small  since  the  rotor  and  wing  contributions  to 
this  derivative  are  of  opposite  signs.  These  considerations  encourage  us  to 
consider  the  assunptlons  meide  for  the  near-hcverlng  regime  to  be  valid  over 
a  considerable  range  at  low  speeds . 

The  trsmsfer  function  for  the  attitude  response  within  this  range  can  be 
expressed  as  follows : 
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(s  -  z„)(s  -  x^)  , _ 

(S  -  Z^)  [B(S  -  Me)(S  -  Xu)  +  gM^ 


(17) 


Hence,  the  mode  Is  cancelled  as  far  as  the  0  response  is  concerned 

(this  mode  is  also  cancelled  for  the  u  resixanse).  This  leads  to  the  follow¬ 
ing  third  order  characteristic  equation  for  low  speeds: 


S(S  -  Mg)(S  -  Xu)  +  gMu  =  0 


(18) 


The  location  of  the  chsuracteristic  roots  can  he  determined  from  the 
following  root  locus  where  the  gain  is  gMu 
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The  larger  gMj^,  the  closer  are  the  oscillatory  roots  to  the  asymptotep 
which  originate  In  1/3(Mq  +  X^).  For  large  enough  the  three  original 

poles  (Mq,  X^,  and  zero)  can  In  effect  be  replaced  by  a  triple  pole  at 
1/3(M§  +  X^).  In  this  case 

“3 


-  1/3(M§  +  X^) 


+  gMu  =  0 


(19) 


can  be  considered  as  an  approximating  characteristic  eq^uatlon.  This  approxl- 
^  matlon  places  the  poles  right  on  the  asymptotes  at  eq.ual  distances  of 

SMu  from  the  1/3(M0  +  X^^)  point.  If  the  actual  roots  are  S3  ,  S3, 
and  If  the  approximating  roots  are  Sf ,  S^,  S*,  the  error  of  approxlmatloh 
can  be  expressed  as  ^ 


(SJ 


S^)  =  2s, 


(S|  -  S2)  =  -  (€^  +  je  ) 


(20) 


(S* 

^3 


S3)  .  . 


Using  the  resiilts  obtained  for  Cy  In  Appendix  II,  It  can  be  shown  that 

1  ^/3(M0  +  Xjj)  -  M§Xu 

-  (21) 


®x«  S 


“\|r 


V 


eM„ 


(22) 


For  a  VTOL  with  a  reasonable  amount  of  static  stability  with  velocity,  M^, 
(gMii'^l/sec“3)  and  expected  values  for  M©  and  X^,  we  find  gec"^. 

Therefore,  for  practical  purposes  It  Is  sufficient  to  consider  only  the 
variations  of  the  parameters  1/3(M§  +  Xu)  and\3  as  affecting  the 

displacement  of  the  roots.  ^ 

=  (Sl)hov  4%  +  Xu)  +  gMu  eJ(60°+l  120°)  ^23) 

1  =  1,  2,  3 

This  means  that  the  roots  will  move  horizontally  to  the  left  for  an  Increas¬ 
ing  damping  and  radially  towards  the  c.g.  of  the  three  roots  for  a  decreas¬ 
ing  static  stability  with  velocity. 
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The  maxlmxjin  speed  to  which  this  approximation  can  he  used  depends  primarily 
upon  the  Increase  of  (=VMy)  with  velocity.  Experimental  results  with 
different  VTOL  planes  (References  1,  2,  3)  show  that  first  Increases  with 
velocity,  because  of  the  inherent  rotor  instability  with  angle  of  attack. 
Shortly  thereafter,  the  tall  stabilizing  contribution  overrides  that  effect 
and  %  decreases  and  becomes  negative.  As  a  result,  there  is  a  velocity 
at  which  the  cubic  characteristic  equation  accurately  represents  the  dy¬ 
namics  of  the  vehicle.  It  is  therefore  reasonable  to  consider  our  approxi¬ 
mation  as  accurate  enovigh  for  the  interval  from  hovering  up  to  this  point. 

2.  NEAR-FORWARD  FLIGHT 


As  the  airplane  approaches  its  forward  flight  configuration,  the  static 
stability  with  angle  of  attack  increases  considerably.  The  transfer 
fimction  for  the  pitch  attitude  becomes 

Is  -  Xu  X^,  0 


-  Zu 

S  - 

0 

-  Mu 

-  Mw 

Ms 

-  Su 

x„ 
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s(s  -  M^) 
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=  Mg 


(S  -  Z^){S  -  Xu)  +  X„Zu 

(S  -  Zj  j^S(S  -  M^)(S  -  Xj  +  gM^J  -  M^,V(S^  -  V  ■  ' 

.  X,S  |z„(s  -  MJ)  -  VM„ 

g/V  ZJ  + 

The  derivative  X^  is  not  negligible  now  in  comparison  with  some  other 
derivatives;  yet  the  X^,  terms  in  both  numerator  and  denominator  can  still  be 
disregarded  on  the  basis  that  at  higher  speeds  Mu  is  small  and  is  much 
larger  than  in  hovering.  In  other  words,  it  is  assvnned 


2o 


Hence 


Vu  Vu 


(s  -  z;^)(s  -  x^) 


(s 


-  Z^)  ^(S  -  M§)(S  -  x^)  +  gM^J  -  %(S^  -  X^S  -  g/v  Z^) 


(25) 


It  Is  a  common  practice  to  study  the  short-period  roots  of  conventional  air¬ 
planes  vising  a  second  order  equation  derived  from  the  fourth  order  charac¬ 
teristic  equation  when  the  higher  frequency  roots  are  far  away  from  the 
other  two  roots.  By  the  same  reasoning,  the  same  procedure  can  he  used  in 
dealing  with  VTOLs  during  the  last  part  of  the  transition;  but  it  should  be 
remembered  that  tha  farther  the  flight  condition  is  from  forward  flight, 
the  poorer  will  be  this  approximation. 

In  order  to  show  the  Influence  of  various  parameters  on  the  approximation, 
the  root  locus  technique  is  used  to  determine  the  characteristic  toots. 
(Reference  4).  First,  the  low-frequency  roots  ("phvigold"  roots)  are 
considered. 

From  the  expression 

S(S  -  Mg)(S  -  X^)  +  gM^  (26) 

which  appears  in  square  brackets  in  the  denominator  of  the  transfer 
fvinction  9/8  (equatlon25)  one  can  determine  three  points  on  the  complex 
plane.  The  roots  of  the  equation  are  on  a  root  locus  as  sketched  in 
Figure  JA  with  the  gain  gM^^. 

Mu  is  assvimed  to  be  small  eind  positive  throughovtt  the  range  near  full  speed. 
Two  of  the  roots  will  then  be  either  on  the  real  axis  within  the  Interval 
(O,  X^)  or  close  to  the  real  axis  on  the  curved  branch,  as  shown  in  the 
sketch.  Next,  another  root  locus  can  be  used  to  obtain  the  roots  of  the 
characteristic  equation.  The  starting  poles  for  this  locus  are  the  three 
poles  obtained  as  roots  of  the  eqviation  (26),  and  a  pole  at  Z^,.  The  ^eros 
are  the  roots  of  the  equation; 

-  X^  -  g/V  Z^  =  0  (27) 

and  the  gain  is  This  locus  is  shown  in  Figure  7B. 
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The  zeros  are  most  likely  a  complex  pair  close  to  the  real  axis .  They  lie 
on  the  vertical  line  passing  throijgh  the  point  (l/2  Gj)  and  their  dis¬ 
tance  from  the  origin  is  (g/v  The  effect  of  a  large  gain,  -  Ma,  is 

to  "bring  the  low-frequency  roots  of  the  characteristic  eqiiation  very  close 
to  the  zeros  given  hy  equation  (27).  The  result  is  the  well-known  approxi¬ 
mation  of  the  phijgold  roots . 

The  above  presentation  shows  how  the  approximation  deteriorates  for  smaller 
In  this  case  the  approximation  overestimates  the  frequency  of  the 
phugold  roots.  Consequently,  (g/v  is  a  limit  fcj'  the  frequency  of 

the  phi^goid  roots.  The  influence  of  all  the  other  derivatives  on  these 
roots  can  be  easily  followed  by  estimating  changes  In  the  geometry  of 
Figure  7. 

The  two  modes  of  motion,  the  ph\agoid  and  the  short -period  mode,  are  well 
separated  during  the  high-speed  phase  of  the  transition.  Since  the  phugoid 
roots  are  considerably  closer  to  the  zeros  of  the  transfer  function  9/6, 
the  short -period  mode  C£in  be  considered  to  be  the  predominant  mode  of  the 
pitch -attitude  response. 

An  analysis  follows  of  the  second  order  approximation  of  the  characteristic 
equation 


-  (M§  +  Z„  +  Xu)>^  +  (Zv  +  Xu)M§  -  +  Z„Xu  =  0  (28) 


For  VTOL  configurations  with  improved  dynamic  and  aerodynamic  characteristics 
the  following  ineq\iality  can  be  expected  to  be  valid  during  the  high-speed 
portion  of  the  transition; 


Me>>  Xu 

This  leads  to  fxirther  simplification  of  the  second  order  equation; 

X^  -  (M^  +  Z„  +  XjX  +  (Z„  +  X^)M§  -1^=0  (29) 


Using  the  root  locxis  technique  is  again  helpful  in  tracing  the  effects  of 
the  derivatives .  The  starting  poles  of  the  root  loctis 
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axe  M§  and  7^  +  3C^.  The  gain  determining  the  roots  is  -  The  roots  lie 

on  a  vertical  line  perpendicular  to  the  real  axis,  crossing  the  axis  at  the 
arithmetic  mean  of  (Zw  +  Xu)  M&. 

The  constant  term  in  eq,viation  (29),  (Zw  +  Xu)M0  -  Hxt  represents  the  square 
of  the  distance  between  either  root  and  the  origin.  If  the  difference  be¬ 
tween  Mg  and  is  considerably  smaller  than  their  absolute  values 

then  their  pr^uct  which  appears  in  the  constant  term  can  be  substituted 
by  the  square  of  their  arithmetic  mean.  (The  error  introduced  by  this  as¬ 
sumption  is  only  6$  if  their  difference  is  half  their  arithmetic  mean  value) 
With  this  siibstlt\itlon  it  can  be  shown  that  (-  %)V2  is  a  good  approxi¬ 
mation  to  the  distance  of  the  roots  from  the  real  axis.  The  real  part  of 
the  short -period  roots  is  approximated,  as  stated  before,  by 
1/2(M0  +  +  X^).  The  dan^jing  ratio  depends  on  the  parameter 

(M^  +  z^  +Xu)yc- 

3.  CQNSTRUCTIOW  OF  TRAHSmCK  LOCI  REAR  HOVERING  AND  NEAR  THE  EMD  OF 
THE  TRANSITION. 


The  approximations  of  the  oscillatory  roots  near  hovering  and  near  the  end 
of  the  transition  can  be  ob^ined  graphically  as  shown  in  Flgrire  8. 
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According  to  the  prevloiis  approximations  two  limped  parameters  in  teyms  of 
stahllity  derivatives  are  introduced  for  the  near-hovering  regime  (  \(  git, 

and  1/3  1^  +  Xy^j  )  and  two  others  for  the  regime  near  the  forward  flight 
end  of  the  transition  (  \j  -  and  l/2  +  2^  +  xj ) .  The  values  of  these 

parameters  can  be  assumed  to  be  known  for  hovering  and  forward  flight 
(completed  transition) . 

The  changes  in  these  parameters  determine  the  approximate  changes  in  the 
oscillatory  roots  near  both  end  points  of  a  transition.  Figure  8  illus¬ 
trates  boundaries  for  assumed  changes  of  the  parameters  with  respect  to 
the  nondimens ional  velocity  fllght>  their  effect  on  varia¬ 

tions  of  the  oscillatory  roots . 

The  assumed  boundaries  are  shown  in  auxiliary  coordinate  systems  by  means 
of  two  lines  for  each  parameter.  Parameter  variations  described  by  any 
line  running  between  these  limit  lines  will  cause  root  variations  within 
the  constructed  resultant  boundaries  on  the  S-plane.  As  an  example,  an 
upper  boundary  line  for  the  location  of  the  short-period  roots  is  obtained 
by  combining  the  line  of  assumed  minima  of  the  damping,  1/2(M0  +  +  X^), 

for  the  near-foiward  flight  situation  with  the  line  of  assumed  maxima  of 
the  freq.uency 

In  the  near-hovering  flight  region  the  roots  are  determined  graphically  by 
their  coordinates  in  an  auxiliary  axis  system  in  which  the  angle  between 
the  axes  is  60°.  The  root  coordinates  in  this  system  eu?e  given  by 
and  1/3 (Mg  +  X^^).  Similarly,  the  upper  and  lower  boundaries,  of  the  roots 
are  determined  ais  above. 

The  plot  of  the  roots  on  the  S-plane  shows  a  gap  between  the  lines  a  a*  and 
b  b '  on  which  the  roots  corresponding  to  0.4  ^  are  located  when  ob¬ 

tained  by  the  near-hovering  and  the  near-forwM^  flight  approximation. 

This  is  the  natural  conseq.uence  of  the  inaccuracies  in  the  approximations. 
The  actual  size  of  the  gap  also  depends  on  the  relative  location  of  the 
hovering  and  the  forweird  flight  roots. 

4.  MID-TRAHSITION  FLIGHT 

The  mid-transition  flight  regime  can  be  defined  as  that  range  where  is 
too  small  for  the  assumption  made  in  the  near-forward-flight  approximation, 
but  cannot  be  neglected  as  was  done  in  the  near-hovering  flight  assumptions. 
Iforeover,  ,  which  depends  upon  the  wing  operating  conditions,  can  change 
largely  within  this  range  whenever  the  wing  operates  in  the  nonlinear  region 
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of  the  lift  cvtrve  vlth  respect  to  effective  angle  of  attack.  Simplifying 
approximations  like  those  made  for  both  end  regions  of  the  transition  re- 
sialt  in  a  grossly  distorted  picture,  because  it  is  characteristic  for  the 
mid-transition  regime  that  all  the  important  stability  derivatives  are  of 
comparable  significance.  Despite  this  fact,  the  construction  of  Flg\are  8 
is  extended  into  the  mid -transition  region.  It  is  interesting  to  note  that 
even  this  oversimplification  gives  an  Indication  of  the  general  behavior  of 
the  oscillatory  roots  thro^lghout  the  transition.  In  partic\i3ar,  both  an 
early  increase  in  and  a  fast  decrease  of  with  increasing  velocity. 
Indicate  a  "dip"  of  the  oscillatory  roots  toward  the  real  axis .  However, 
it  cannot  be  predicted  by  the  illustrated  procedure  how  far  this  dip  actu¬ 
ally  extends  in  mid -transit ion  flight.  In  fact,  two  alternatives  for  the 
general  behavior  of  the  roots  can  be  recognized.  It  is  possible  that  the 
hovering  roots  move  (perhaps  throijgh  a  "dip")  toward  the  short-period 
roots  and  that  the  two  real  roots  in  hovering  combine  into  the  phugold 
roots  in  forward  flight  (Figure  9A).  The  alternative  possibility  is  that 
the  hovering  roots  move  toward  the  phugold  roots  during  the  transition  and 
that  the  real  roots  in  hovering  combine  and  break  away  toward  the  short- 
period  roots  somewhare  in  mid-transition  (Figure  9B)»  The  loci  in 
Figure  9  are  not  root  loci  in  the  \isual  sense  but  describe  the  motion  of 
the  characteristic  roots  in  the  congslex  plane  during  transition. 

Further  application  of  the  root  locvis  techniq,ue  indicates  that  these  two 
different  behaviors  for  the  roots  depend  primarily  on  the  relationship  be¬ 
tween  Zw  and  during  mid-transition.  If  Met  increases  much  earlier  in  the 
cotirse  of  the  transition  than  Zw,  the  loctis  of  Flg\are  9A  will  be  the  re¬ 
sult.  If  Zw  Increases  earlier  than  Vcc,  the  locus  will  be  of  the  type 
ill\istrated  in  Figure  9®* 

In  general,  the  variation  of  ,  Mg  ,  Z^  ,  and  in  mid-transition  de¬ 
termines  the  eventual  amovuat  of  "dip"  in  the  movement  of  the  oscillatory 
roots.  A  sharp  dip  Implies  a  large  decrease  and  Increase  within  a  rela¬ 
tively  small  speed  range.  In  such  a  case  the  resiilts  of  the  previous 
section  on  the  effects  of  fast  time  variation  apply. 

Apparent  deviations  from  the  frozen  system  approximation  can  be  expected 
to  Influence  the  handling  qualities  and  feedback  control.  It  has  been  shown 
that  a  fast  decrease  in  frequency  near  the  Imaginary  axis  can  ca\ase  a  con¬ 
siderable  decrease  in  apparent  damping.  The  approximate  analysis  presented 
It.  this  section  has  identified  the  possible  causes  leading  to  this  situation. 
The  method  can  be  used  in  a  quick  evaluation  of  the  effects  of  design 
changes  on  the  behavior  of  the  oscillatory  characteristic  roots  during  a 
transition. 
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5.  OBCILLATOBY  ROCCTS  IN  TRANSITION 


In  order  to  Illustrate  transitions  of  the  simplified  second  order  model 
described  earlier,  the  piecewise  straight  transition  loc\is  shown  in 
Figxire  IQA.  has  been  chosen.  The  component  of  speed  of  variation  parallel 
to  the  real  axis  was  kept  consteuat.  Pulse  inputs  were  applied  automatically 
at  the  instants  when  the  roots  were  in  the  positions  marked  1,  2,  3# 
only  one  point  was  tested  in  each  transition.  Transients  for  the  fairly 
fast  total  transition  time  of  32  seconds  are  illustrated  in  Figures  ICB 
through  ICE.  This  time  is  approximately  five  times  the  period  of  the  hover¬ 
ing  or  mid-transition  response.  The  frozen  system  piilse  response  (a)  and 
the  responses  for  the  accelerating  (b)  and  decelerating  (c)  transition 
are  shown  for  comparison  in  each  figure.  Both  the  displacement  and  the 
velocity  responses  are  shown. 

Transients  for  an  even  faster  transition  of  l6  seconds  (approximately  2.5 
periods)  are  shown  superimposed  in  Figures  UA  to  ITD  for  the  points  illus¬ 
trated  in  Figure  IQA.,  Note  the  considerable  destabilizing  effect  during 
deceleration  and  the  stabilizing  effect  during  acceleration. 
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FEEDBACK  CQNSIDEBATIQNS 


In  this  final  section  of  the  report  some  aspects  of  artificial  pltch- 
•^ttitude  stability  a\agmentation  of  VTOL  aircraft  are  discussed.  These 
aspects  are  quite  general  in  their  nature  and  part  of  the  discussion  ap¬ 
plies  to  all  time  varying  feedback  control  systems  In  general.  Some 
important  problems  are  omitted,  such  as  the  provision  of  sufficient  con¬ 
trol  authority  for  the  feedback  system.  Emphasis  Is  placed  on  those 
aspects  which  are  connected  with  the  time  varying  natiire  of  the  problem. 

The  range  of  variations  can  be  estimated,  based  on  the  frozen  approxi¬ 
mation.  The  deviation  from  the  frozen  can  be  estimated,  based  on  the 
resiilts  presented  previously  in  this  report. 

In  order  to  provide  a  basis  for  the  choice  of  the  feedback  configurations, 
the  goal  to  be  achieved  miist  first  be  determined.  Since  the  concern  Is 
with  piloted  aircraft,  some  variation  in  frequency  as  well  as  in  damping 
of  the  predominant  closed  loop  response  can  be  tolerated  because  of  pilot 
adaptivity.  A  range  of  low  frequencies  within  a  2;1  ratio  and  a  range  of 
relative  damping  ratios  between  .2  and  :,7  can  be  handled  easily  by  a  pilot 
Such  tolerances  define  a  considerable  area  in  the  complex  plane  within 
which  the  predominant  closed  loop  roots  may  possibly  be  located. 

We  assume  a  basic  rate  and  attitude  feedback  loop  in  which  the  feedback 
gains  have  to  be  adjusted  because  of  the  time  variation  of  the  system. 

The  following  alternatives  will  be  considered  as  fundamentally  differing 
feedback  adjustments  for  time  varying  feedback  systems . 

1.  ADJUSTMEMT  BY  DIRECT  FEEDBACK 

In  principle,  it  would  seem  desirable  to'make  adjustments  in  the  basic  feed¬ 
back  loop  according  to  the  variation  of  the  variable  which  is  the  physical 
cause  of  the  time  variation.  This  is  the  dynamic  pressure  or  velocity  for 
VTOLs  and  helicopters,  and  a  more  complicated  pressure  distribution  in  the 
case  of  ground  effect  machines.  If  the  time  varying  system  is  well  known, 
then  the  proper  gain  adjustments  as  fvinctlons  of  the  velocity  can  be  pre¬ 
determined.  The  instrumentation  of  such  a  scheme  involves  the  measurement 
of  the  velocity,  some  nonlinear  device  determining  the  amount  of  gain  adjust¬ 
ment  and  the  gain  adjusting  element  itself  (Figure  12A).  Several  dlffl- 
cixltles  make  this  approach  unfeasible  for  practical  application  to  low- 
speed  flying  machines.  Dynamic  pressures  cannot  be  meas\ared  well  enough  at 
very  low  speeds .  However,  even  if  they  could  be  determined  adequately,  two 
other  fundamental  problems  remain.  In  the  frozen  system  approach,  the  un¬ 
certainty  about  the  actual  variations  of  derivatives  makes  the  determination 
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of  the  nonlinear  adjustment  functions  saiffiwhat  arbitrary.  Additional  con¬ 
ditions  such  as  the  amoiint  of  acceleration  or  deceleration  at  a  certain 
velocity  cause  large  variations  and  additional  comlexlty. 

The  adjustment  by  direct  feedback  means  an  additional  feedback  loop  around 
the  system  with  the  basic  loop.  The  two  loops  are  coupled  through  the  gain 
adjiostment.  In  the  case  of  slow  time  variation,  the  variable  on  which  the 
adjustment  depends,  the  velocity,  can  be  filtered  in  order  to  extract  the 
average  variation  from  the  time  history  of  the  velocity.  For  fast  time 
variations,  however,  when  significant  changes  take  place  within  one  period 
of  a  predominant  mode,  the  filtering  time  constant  has  to  be  so  short  that 
dynamic  coupling  between  the  two  loops  cannot  be  avoided.  This  causes  an 
additional  problem  because  the  dynamics  of  the  entire  system  become  in¬ 
creasingly  complicated.  Actually,  the  linearization  which  leads  to  a  time 
varying  linear  system  is  not  valid  in  this  case.  For  these  reasons  gain 
adjustments  by  means  of  direct  feedback  should  not  be  used  for  the  purpose 
in  mind  here. 

2 .  ADAPTIVE  FEEDBACK 

The  principle  of  adaptive  feedback  was  Introduced  to  handle  control  pro¬ 
blems  where  the  system  to  be  controlled  is  not  sufficiently  well  defined 
for  a  conventional  feedback  system  synthesis.  The  term  "adaptive”  is  used 
here  for  systems  in  which  the  performance  of  the  feedback  control  is  de¬ 
termined  in  terms  of  a  performance  index,  and  adjiistments  are  made  accord¬ 
ingly  in  the  basic  feedback  loop  (Figure  12B).  Based  on  the  frozen  sys¬ 
tem  approximation,  two  major  difficulties  which  were  connected  with  the 
first  approach  could  be  eliminated  by  an  adaptive  control  system  but  at  the 
expense  of  added  complexity.  These  R.rc  the  problems  of  velocity  measvure- 
ment  at  low  speeds  eind  the  effects  of  design  uncertainties.  However,  the 
less  the  time  available  for  system  identification,  the  larger  will  be  the 
instantaneovis  deviations  of  the  gain  settings  from  the  optimum.  If  an 
adaptive  system  is  expected  to  follow  changes  within  one  period  of  a  pre¬ 
dominant  mode,  then  dynamic  coupling  between  the  adaptive  loop  ajid  the 
basic  feedback  loop  cannot  be  avoided.  Fast  system  identification  is  com¬ 
plicated  even  theoretically  and  could  lead  to  very  complex  instr\imentatlon. 

One  way  of  getting  around  this  problem  is  to  suppress  the  original  pre¬ 
dominant  mode  of  the  system.  This  Implies  that  a  considerably  higher  feed¬ 
back  gain  can  be  applied  than  would  be  necessary  for  stability  avigmentation 
itself,  so  that  the  natural  freq.uency  of  the  predominant  closed  loop  mode 
is  Increased  to  be  fast  with  respect  to  the  time  variation.  The  resvilt 
is  a  very  tight  control  system  implying .  that  more  control  authority  must 
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be  available  than  woiild  be  necessary  for  stability  augmentation  only.  For 
example,  in  the  Minneapolis -Honeywell  adaptive  autopilot  such  a  tight  loop 
arovuad  the  aircraft  is  used  to  make  the  aircraft  follow  a  dynamic  model 
response . 

Whenever  the  application  of  adaptive  control  is  considered,  as  much  a  priori 
knowledge  as  possible  should  be  used  in  order  to  make  the  design  as  slinple 
and  as  efficient  as  possible.  It  has  been  shown  in  the  section,  "Approxi¬ 
mations  to  the  Longitudinal  Dynamics  of  a  VTOL  Aircraft,"  that  even  a  very 
simplified  description  leads  to  considerable  information  on  the  behavior  of 
the  predominant  roots .  Determining  the  roots  in  the  forward  flight  configu¬ 
ration  and  hovering  are  classical  problems .  An  area  inside  which  the  oscil¬ 
latory  roots  move  in  transition  can  be  reasonably  well  approximated;  the 
actual  root  will  not  be  too  far  from  the  crude  approximation.  A  q.\aestion 
must  be  raised  as  to  whether  the  complexity  of  an  eidaptive  system  is  Justi¬ 
fied  or  whether  a  much  simpler  approach  might  not  be  as  effective. 

3.  PROGRAMMEID  ADJTJSTMEaCTS 

The  considerable  amount  of  a  priori  knowledge  and  the  rather  broad  toler¬ 
ance  permitted  becai:ise  of  the  pilot’s  adaptivity  stiggest  the  feasibility 
of  simple  feedback  gain  adjustments.  The  range  of  open  loop  damping  is 
determined  by  the  hovering  and  the  fojrward  flight  roots.  The  frequency 
range  is  influenced  largely  by  the  "dip"  of  the  transition  locus  (Figure  8). 
The  simplifying  approximations  made  throughout  this  report  can  serve  as  a 
basis  for  programming  the  feedback  adjustments.  Programming  means  that  the 
adjustments  are  "open  loop"  adjustments  rather  than  being  derived  from  the 
output  of  the  system  as  in  the  previous  two  configurations  (Fig\ire  12C). 
Adjustments  are  programmed  as  functions  of  the  wing  angle  in  the  Tri- 
Service  VTOL  aircraft.  The  closed  loop- roots  of  an  approximating  model 
can  be  confined  to  a  small  region  of  the  complex  plane  by  means  of  rather 
simple  gain  programs.  The  neglected  poles  and  zeros,  variations  of  the 
c.g.  of  the  aircraft,  and  deviations  due  to  acceleration  and  deceleration 
will  cause  the  closed  loop  roots  to  move  within  a  larger  region  than  the 
one  predicted  by  neglecting  these  effects .  Thus  it  is  necessary  to  rely  on 
the  adaptivity  of  the  pilot.  The  effect  of  fast  time  variation  must  also 
be  taken  into  account,  using  the  relationships  established  in  the  Section, 
"Analog  Coii5)uter  Study"  and  Appendix  I.  It  has  been  shown  that  the  apparent 
effects  due  to  fast  time  variation  are  strongest  near  the  imaginary  axis. 
Therefore,  with  increased  closed  loop  damping  the  effect  of  the  fast  time 
variation  becomes  less  api)arent. 
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The  simplicity  of  programmed  feedback  adjvistments  Is  very  attractive. 
Further  Investigation  of  the  sensitivity  of  such  a  system  to  variations 
and  uncertainties  of  system  parameters  is  imder  way  in  order  to  establish 
the  full  potential  and  limitations  of  this  approach.  Whenever  the  per¬ 
formance  of  this  approach  can  be  made  acceptable  it  should  be  preferred 
to  the  more  elegant  b\jt  much  more  complex  adaptive  approach. 

Similar  conclusions  were  reached  in  a  previous  report  (Reference  5)  in 
connection  with  the  artificial  stability  augmentation  of  ground  effect 
machines . 
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APPENDIX  I 


NCRMALIZATIQN 


The  normalization  of 

pt^t)x  +  +  ^to^t)  =  0  (Al) 

through  y  =  and  t'  =  Wpt  where  ojp  =  u>p^ 

leads  to 

y  +  Pt'y  +  (1  +  QLt')y  =  0  (A2) 

This  transformation  is  eqxiivalent  to  a  horizontal  shift  of  the  frozen  poles 
by  Imaginary  axis  and  a  scale  expansion  on  both  axes  by  the 

factor  l/up.  Lines  passing  through  the  frozen  point  on  the  plane  eure 
transformed  into  lines  passing  through  the  (+j,  -j)  point  on  the  Sy  plane. 
The  transformation  is  angle  preserving.  The  linear  time  variation  of  the 
coefficients  causes  the  frozen  roots  to  move  along  circles  centered  on  the 
real  axis.  These  circles  go  over  into  circles  in  the  course  of  the  trans¬ 
formation.  The  relationships  between  P,  Q  and  p,  p.  was  given,  by  definition 

P=P;  Q  =  p  +  p  tan  T|f  (>k  =  tan"^  )  (A3) 

-  4 

The  slope  of  the  variation  as  well  as  the  rate  of  variation  in  the  normall- 
zed:‘5y  plane  can  be  determined  as  functions  of  p  and  p.  The  normalized 
frozen  locus  is  obtained  by  taking  t'  as  a  parameter  in  the  following 
characteristic  eq\iatlon. 


X  +  (2Cp‘^o  + 


+  Pt'Sy  +  (1  +  Qt')  =  0 


Sy(t')  =  1/2  j-  Pt'  ±  \ 


r 


(Pt')^  -  4(1  +  Qt') 


For  the  variation  of  Sy  in  the  vicinity  of  t’  =  0 


<iSy(f) 

"IP 


=  l/2(P  ±  JQ) 

fsO 


(A4) 

(A5) 


(  A6) 
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The  slope  ot  the  locvis  Is  detennlned  by 


dSy(t'  ) 

"dP 


t'=0 


=  7  =  tan  i  =  tan“^  - 


P 


\F 


(AT) 


=  tan  ^  where  \|f  Is  the  angle  enclosed  by  the  radius  vector 


\l 


1  -  it 


of  the  pole  and  the  imaginary  axis,  p/p  can  be  considered  the  tangent  of 
a  fictitious  angle  ♦  .  This  permits  the  Illustration  of  the  full  range 
of  p  and  p  ,  including  p  =  0,  on  a  finite  scale.  The  angle  7  of  the 
normalized  frozen  root  variation  is  related  to  p/p  and  of  the  original 
root  through  the  following  relationship 


tan  7  =  tan  ♦+  tan  where 


♦  =  tan”^  Z. 


(AS) 


This  relationship  is  plotted  with  as  a  parameter  on  Figure  13.  The 

speed  of  root  variation  is  also  determined  by  p,  p  and 

dSy(f) 

"dp 


t'=0 


»  i/pS^  P^  + 


“TP - 

V  -  Ci 


(A9) 


For  poles  along  the  imaginary  axis,  where  =  0,  this  becomes 

1/2  p^  +  p2 

The  general  case  can  be  related  to  this  value. 

'aSy(t')| 


dt' 


t'=o 


=  l/2\p=  ^  \ 


P)" 


=  1/2 


.\l 


N 


a2  2 

P  +  >1 


(AlO) 


a2  2 

P  +  P 


The  coefficient  of 
in  Figure  14- . 


l/2\ 


\jl  +  (tan  ♦+  tan  'ffY 
^  1  +  (tan  ♦ )^ 


p 

P  +  p  with  as  a  parameter,  is  Illustrated 
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The  shapes  of  the  transformed  frozen  pole  loci  with  time  variation  are 
circles  with  their  centers  on  the  real  axis  at  the  distance 


Q/P 


tan  '!>  +  tan  \lr 


« 


3h 


APPEMDIX  II 


NEAR-HOVERING -FUGHT  APPROXIMATION 

We  can  determine  the  error  made  In  the  determination  of  roots  when  the  near- 
hovering  characteristic  eq.mtlon 

S(S  -  M§)(S  -  X^)  +  gM^  =  0  (All) 

is  substituted  by  the  following  equation 

1^  -  l/3(Me  +  x^)  ^  +  gM^  =  0  (A12) 

For  convenience  we  shift  the  origin  of  the  coordinate  system  into  the 
c.g.  of  the  roots  of  both  equations  into  the  (1/3K9  +  0)  point.  This 

transformation  does  not  change  the  errors  nor  the  c.g.  of  the  roots. 

The  new  cubic  equation  obtained  from  equation  (All)  is  of  the  type 

+  BX  +  c3  =  0  (AI3) 
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It  can  also  te  shown  that  B  <  C 


We  solve  the  approximating  equation  which  is  obtained  from  equation  (A12) 
as 


(Al4) 


then  add  to  the  roots  correcting  terms  which  are  explicit  functions  of  the 
coefficients  B  and  C.  The  roots  of  equation  (Al4)  are 

^1  =  C,  4^3  =  1/2C  +  Jl/2  \\3  C  (A15) 

These  roots  form  an  equilateral  triangle  in  the  complex  plane.  The  roots 
of  equation  (A13)  differ  from  these  roots.  The  components  and  Sy  of 

the  error  are  indicated  in  the  following  figiire 
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The  roots  of  the  exact  eq.\iation  can  he  expressed  as 


Xi  =  -  C  -  26^  ,  ^2^3  =  1/2C  +  ±j(l/2  -\|C  -  €^)  (ai6) 

The  coefficient  B  is  eq.ijal  to  the  sum  of  all  double  products  of  the  roots 


B  =  (-  C  -  2e^)2  (l/2C  +  e^)  -  — 


C  -  2e, 


(The -product  of  the  two  Imaginary  roots  is  expressed  in  terms  of  the 
negative  triple  product  of  the  roots  -  c’  divided  by  the  real  roots . ) 

Using  series  expansion,  and  neglecting  the  terms  which  are  higher  than 
second  order  in  2ex/C 


p  OP  2€y  ^€x 

B  «  -  C  -  l^«^c  -  ke^+  C^(l  -  ) 


(A17) 


‘  'Sc 

Similarly,  the  coefficient  can  be  expressed  as  the  product  of  the  roots 

as  follows 


-  C3  =  (.  C  -  2€x) 


(1/2C  + 


€x)  +  (1/2  -3  C  -  B  f 


(A19) 


+  c2(  -  3ex  +  ^  3|~) 


In  terms  of  the  stability  derivatives  we  have 


(A20) 


Cy  i  (l/3(M0  +  -  MeX^ 
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Finally 


arc  tan  ^  r:?  60° 


(A21) 


These  results  indicate  that  the  complex  roots  of  the  exact  eq.uatlon  (A13) 
lie  on  lines  inclined  60^^  with  respect  to  the  horizontal  and  pass  through 
the  complex  roots  of  the  approximating  eq^uation  (Al4).  The  distance  be¬ 


tween  each  of  these  roots  and  the  related  roots  of  the  exact  equation  (A13) 
is  approximately  equal  to  2 


These  conclusions  yield  a  new  and  quick  method  for  estimating  the  roots  of 
a  cubic  equation  of  the  type  given  by  (A13). 
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APPENDIX  III 


votrs 


FIGURE  2.  TYPICAL  ISE  SELECTIVITY 


FIGURE  3.  APPROXIMATING  ROOTS 
SPEED  OF  VARIATION  r^O*  +P*-^ 
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FIGURE  4.  approximating  ROOTS 


FIGURE  5.  NORMALIZED  APPROXIMATING  ROOTS 


FIGURE  6.  BOUNDARY  BETWEEN  INCREASED 
AND  DECREASED  APPARENT  DAMPING 


FIGURE  9.  TWO  FOSSISLE  ROOT  TRANSITIONS 


FICUPTE  II.  !•  SECOND  TRANSITION 


FieUflE  MD.  TRANSIENTS  IN  POINT  4 


-80 

r  IN 

DEGREES 

^  IN 

DEGREES 

*  tan"' 


y  *  tan"' 


iJL-Sl— 

^  4^ 


) 


FIGURE  13.  SLOPE  OF  FROZEN  LOCUS 


FIGURE  14. 
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variable  feedback  configurations:  direct  feedback  variable  feedback  configurations:  direct  feedback 

adjustments,  adaptive  feedback,  and  programmed  feedback  adjustments,  adaptive  feedback,  and  programmed  feedback 
adjustments .  adjustments . 
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